Recent astrophysical and accelerator based results on the Hadronic Equation of State 
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In astrophysics as well as in hadron physics progress has recently been made on the determination 
of the hadronic equation of state (EOS) of compressed matter. The results are contradictory, 
however. Simulations of heavy ion reactions are now sufficiently robust to predict the stiffness of 
the (EOS) from (i) the energy dependence of the ratio of K + from Au+Au and C+C collisions 
and (ii) the centrality dependence of the K + multiplicities. The data are best described with a 
compressibility coefficient at normal nuclear matter density k around 200 MeV, a value which is 
usually called "soft" The recent observation of a neutron star with a mass of twice the solar mass 
is only compatible with theoretical predictions if the EOS is stiff. We review the present situation. 
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How much energy is needed to compress nuclear matter? The answer to this question, the determination of 
E/A(p,T), the energy /nucleon in nuclear matter in thermal equilibrium as a function of the density p and the 
temperature T, has been considered since many years as one of the most important challenges in nuclear physics. 
This quest has been dubbed "search for the nuclear equation of state (EoS)". 

Only at equilibrium density, po, the energy per nucleon E/A(p = po, T = 0) = — 16 MeV is known by extrapolating 
the Weizsacker mass formula to infinite matter. Standard ab initio many body calculations do not allow for a 
determination of E/A(p, T) at energies well above the saturation density because the low density many body expansion 
schema ( Bruckner G- matrix) breaks down and therefore the number of contributing terms is exploding. Therefore 
in nuclear reaction physics another strategy has been developed. Theory has identified experimental observables in 
nuclear reaction physics or in astrophysics which are sensitive to E/A(p,T). Unfortunately these observables depend 
as well on other quantit ies which are either unknown or little known (like cross sections wit h resonance in the entrance 
channel) or difficult to isses theoreticalljy (like the resonance lifetinles in hot and dense matter). It was hoped that 
comparing many observ; iblesJbr different systems and different energies with the theoretic; 1 predictions these unknown 
or little known quantities c£,n be determined experimentally and that finally the depenience of the observables on 
E/A(p,T) can be isolated. KA([/ + ) = 

In astrophysics the m .clea : EoS plays an implrtarklmle in bmary /trargers involving h ack holes and neutron stars 
[1], in double pulsar [2 , in the mass-radius relation of neutron stars [3, 4] and in supe -novae explosions [5]. For a 
recent review on {fy© t< spies we refer to [6] . Unfortunately, as in nuclear reaction physics; there are always other little 
known processes or proiertSJs which have to be understood before the nuclear EoS dependence can be isolated. We 
discuss here as example of t ie mass-rad.us relation of neutron stars. Fig. 1 shows the i eutron star masses in units 
of the solar mass for di iierejit types of binaries. These masses axe concentrated at aroind 1-1.5 solar masses. Fig. 
2 shows a theoretical p -edictRni of the mass-radius relation for i eutron stars using different EoS. Since the nature 
of the interior of neutrcn stars lj^iot krown (in contradiction to what the name sugges ;s) one may suppose that it 
consisSSf hadrolis dr o ^quarks. But^vjin if it §]Qr|t^of hadron s there are speculations that there is a K~ or a 7r~ 
condensate or thai? tnen are hyperons i3(k«ujilibrium with nuclear re sonances. The same i 3 true if the interior consists 

of quarks. Little known color-flavor I I I qua*W modify the EoS at densities which are reached in the 

interior of the neutron star. For a detailed discussionoi*SS^jes ; phenomena we refer o ref. [6]. We see that the 
observed masses of neutron stars are ccOaatible with almost all <p''lh4^r hadron based" 3oS as long as the radius is 
unknown. Radii, howev 3r, are very difficultto *jj£asure. Because 3iinilarprWi4mis appea r also for other observables, 
up to recently the astro ahysical observations of neft»ftn stars did not help mucrira0arr' )w down the uncertainty on 
the nuclear EoS.^ 2 — — 

This situation lias ch anged dramatically in the last year with ffie oTSServcTBbrHlP a r eutron star with a mass of 
two solar masses [7]. I ' this observation is finally confirmed the mass/radius predictioi of fig. 2 excludes that the 
interior of a neutron sta As made by quarks [4], (sY^r^!^ nuclear EoS, which will be defiled below, will be excluded. 
This is confirmed by tf e calculation of Maieron [8] which uses a MIT bag model or a color dielectric models EoS 
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FIG. 1: Measured and estimated masses of neutron stars in radio binary pulsars and in x-ray accreting binaries. Error bars are 
la. Vertical dotted lines show average masses of each group (1.62 M@, 1.34 Mq and 1.56 Mq); dashed vertical lines indicate 
inverse error weighted average masses (1.48 Mq, 1.41 M Q and 1.34 Mq). The figure is taken from ref [4] 



to describe the quark phase. Baldo [9] argue that this conclusion may be premature because it depends too much 
on the equation of state of the quark phase. If one replaces the MIT bag model equation of state by that of the 
Nambu - Jona-Lasinio (NJL) Lagrangian under certain conditions (no color conducting phase) larger masses may be 
obtained. The standard NJL Lagrangian lacks, however, repulsion and in view of the momentum cut-off, necessary 
to regularize the loop integrals, and the coupling constants in the diquark sector, which are not uniquely determined 
by the Fierz transformation, quantitative prediction at high quark densities are difficult in this approach even if 
qualitative agreement with pQCD calculation can be found [10]. 

In heavy ion reactions three observables have been identified which are - according to theoretical calculations - 
sensitive to E/A(p,T) at densities larger than p : (i) the strength distribution of giant isoscalar monopole resonances 
[11, 12], (ii) the in-plane sidewards flow of nucleons in semi-central heavy ion reaction at energies between 100 A 
MeV and 400 A MeV [13] and (iii) the production of K + mesons in heavy ion reactions at energies around 1 A GeV 
[14]. We will discuss these approaches below. Although theory has predicted these effects qualitatively, a quantitative 
approach is confronted with two challenges: a) The nucleus is finite and surface effects are not negligible, even for the 
largest nuclei and b) in heavy ion reactions the reacting system does not come into equilibrium. Therefore complicated 
non-equilibrium transport theories have to be employed and the conclusion on the nuclear EoS can only be indirect. 

(i) The study of monopole vibrations has been very successful, but the variation in density is minute. Therefore, 
giant monopole resonances are sensitive to the energy which is necessary to change the density of a cold nucleus 
close to the equilibrium point po- According to theory the vibration frequency depends directly on the force which 
counteracts to any deviation from the equilibrium and therefore to the potential energy. The careful analysis of the 
isoscalar monopole strength in non-relativistic [11] and relativistic mean field models has recently converged [12] due 
to a new parametrization of the relativistic potential. These calculations allow now for the determination of the 

compressibility k = 9p 2 d E { J p p P,T ^ \p=p which measures the curvature of E/A(p,T) at the equilibrium point. The 
values found are around k = 240 MeV and therefore close to what has been dubbed "soft EoS" . It agrees as well 
with the prediction of nuclear matter calculations based on nucleon-nucleon scattering data [15] which give a value of 
about k = 250 MeV. 

(ii) If the overlap zone of projectile and target becomes considerably compressed in semi-central heavy-ion collisions, 
an in-plane flow is created due to the transverse pressure on the baryons outside of the interaction region with this 
flow being proportional to the transverse pressure. In order to obtain a noticeable compression, the beam energy has 
to be large as compared to the Fermi energy of the nucleons inside the nuclei and hence a beam energy of at least 
100 A MeV is necessary. Compression goes along with excitation and therefore the compressional energy of excited 
nuclear matter is encoded in the in-plane flow. It has recently been demonstrated [16] that transport theories do not 




FIG. 2: Mass-radius diagram for neutron stars. Black (green) curves are for normal matter (SQM) EoS [for definitions of the 
labels, see [4]]. Regions excluded by general relativity (GR), causality and rotation constraints are indicated. Contours of 
radiation radii Roo are given by the orange curves. The figure is from [4]. 



agree quantitatively yet and therefore former conclusions [17] have to be considered as premature. 

(iii) The third method is most promising for the study of nuclear matter at high densities [18] and I will discuss 
it in detail this talk. K + mesons produced far below the NN threshold cannot be created in first-chance collisions 
between projectile and target nucleons. They do not provide sufficient energy even if one includes the Fermi motion. 
The effective energy for the production of a K + meson in the TV TV center of mass system is 671 MeV as in addition 
to the mass of the kaon a nucleon has to be converted into a A to conserve strangeness. Before nucleons can create a 
K + at these subthreshold energies, they have to accumulate energy. The most effective way to do this is to convert 
first a nucleon into a A and to produce in a subsequent collision a K + meson via AJV — > NK + A. Two effects link the 
yield of produced K + with the density reached in the collision and the stiffness of the EoS. If less energy is needed 
to compress matter (i) more energy is available for the K + production and (ii) the density which can be reached 
in these reactions will be higher. Higher density means a smaller mean free path and therefore the A will interact 
more often increasing the probability to produce a K + and hence it has a lower chance to decay before it interacts. 
Consequently the K + yield depends on the compressional energy. At beam energies around 1 A GeV matter becomes 
highly excited and mesons are formed. Therefore this process tests highly excited hadronic matter. At beam energies 
> 2 A GeV first-chance collisions dominate and this sensitivity is lost. 

For the third approach different transport theories have converged. This was possible due a special workshop at 
the ECT* in Trento/Italy where the authors of the different codes have discussed their approaches in detail and 
common solutions have been advanced. The results of this common effort have been published in [19]. As an example 
we display here the K + p t spectra at midrapidity obtained in the different transport theories at different energies. 
Because with each K + N rescattering collision the slope of the K + spectra changes the slope of the pt spectra encodes 
not only the K + momentum distribution at the time point of production but also the distribution of the number of 
rescatterings. It is therefore all but trivial. Without the KN potential the slopes are almost identical and even the 
absolute yield which depends on a correct modeling of the Fermi motion of the nucleons is very similar. If we include 
the KN interaction which is not identical in the different approaches (see [19]) we still observe a very similar slope 
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PIG. 3: Final K + transverse momentum distribution at 6=1 fm, \y cm < 0.5 1 and with an enforced A lifetime of 1/120 MeV 
(top row without, bottom row with KN potential) in the different approaches [19]. 

for most of the programs. 

Due to this progress the simulation programs can now be used to extract yet theoretically inaccessible information 
like the hadronic EOS [18]. Two independent experimental observables, the ratio of the excitation functions of the 
K + production for Au+Au and for C+C [20, 21], and a new observable, the dependence on the number of participants 
of the K + yield show that nucleons interact with a potential which corresponds to a compressibility of k < 200 MeV 
in infinite matter in thermal equilibrium. This value extracted for hadronic matter at densities around 2.5 times the 
normal nuclear matter density is very similar to that extracted at normal nuclear matter density. A key point here 
is to demonstrate that the different implementation of yet unsolved physical questions, like the NA — > K + AN cross 
section, the KN interaction as well as the life time of the nuclear resonances in the hadronic environment do not 
affect this conclusion. 

In order to determine the energy which is necessary to compress infinite nuclear matter in thermal equilibrium by 
heavy ion reactions in which no equilibrium is obtained one chooses the following strategy: The transport theory 
calculates the time evolution of the quantal particles described by Gaussian wave functions. The time evolution is 
given by a variational principle and the equations one obtains for this choice of the wave function are identical to 
the classical Hamilton equations where the classical two-body potential is replaced by the expectation value of a 
3 parameter Skyrme potential. The Skyrme potential is a simple approximation to the real part of the Bruckner 
G-matrix which is too complicated for performing simulations of heavy ion collisions. For this potential the potential 
energy in infinite nuclear matter is calculated. To determine the nuclear EoS we average this (momentum-dependent) 
two-body potential over the momentum distribution of a given temperature T and add to it the kinetic energy. 
Expressed as a function of the density we obtain the desired nuclear EoS E /A(p,T). Our two-body potential has 
five parameters which are fixed by the binding energy of infinite nuclear matter at po , the compressibility k and the 
optical potential which has been measured in pA reactions [22]. 

Once the parameters are fixed we use the two-body potential with these parameters in the transport calculation. 
There is an infinite number of two-body potentials which give the same EoS because the range of the potential does 
not play a role in infinite matter. The nuclear surface measured in electron scattering on nuclei fixes the range, 
however, quite well. The uncertainty which remains is of little relevance here (in contradiction to the calculation of 
the in-plane flow which is very sensitive to the exact surface properties of the nuclei and hence to the range of the 
potential). 

We employ the Isospin Quantum Molecular Dynamics (IQMD) with momentum dependent forces. All details of the 
standard version of the program may be found in [22]. The standard version is supplemented for this calculation with 
all inelastic cross sections which are relevant for the K + production. For details of these cross sections we refer to 
[23] . Unless specified differently, the change of the K + mass due to the kaon-nucleon (KN) interaction according to 
m K (p) = ttIq(1 — 0.075^) is taken into account, in agreement with recent self-consistent calculations of the spectral 
function of the K + [24]. The A potential is 2/3 of the nucleon potential, assuming that the s quark is inert. The 
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calculations reproduce the experimental data quite well as can be seen in fig. 4 where we compare the experimental 
and theoretical K + spectra for different centrality bins for 1.48 AGeV Au+Au. This figure shows as well the influence 
of the K + N potential which modifies not only the overall multiplicity of K + due to the increase of the in medium mass 
but also the spectral form confirming the complexity of the transverse momentum spectrum. In order to minimize 
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FIG. 4: K + spectra for different centrality bins as compared with (preliminary] experimental data from the KaoS collaboration 

the experimental systematical errors and the consequences of theoretical uncertainties the KaoS collaboration has 
proposed to study not directly the excitation function of the K + yield but that of the yield ratio of heavy to light 
systems [20]. Calculations have shown that ratios are much less sensitive to little known input parameters because 
they affect both systems in a rather similar way. We have shown in fig. 4 that the absolute yields are well reproduced 
in our simulations. Therefore we can use this ratio directly for a quantitative comparison with data. These ratios are, 
on the contrary, quite sensitive to the nuclear potentials because the compression obtained in the Au+Au collisions 
is considerable (up to 3po) and depends on the nuclear EoS whereas in C+C collisions the compression is small and 
almost independent on the stiffness of the EoS. 

Figure 5 shows the comparison of the measured ratio of the K + multiplicities obtained in Au+Au and C+C reactions 
[20] together with transport model calculations as a function of the beam energy. We see, first of all in the top row, 
that the excitation function of the yield ratio depends on the potential parameters (hard EoS: k = 380 MeV, thin 
lines and solid symbols, soft EoS: k = 200 MeV , thick lines and open symbols) in a quite sensible way and - even 
more essential - that the prediction in the standard version of the simulation (squares) for a soft and a hard EoS 
potential differ much more than the experimental uncertainties. The calculation of Fuchs et al. [21] given in the same 
graph, agrees well with our findings. 

This observation is, as said, not sufficient to determine the potential parameters uniquely because in these transport 
theories several not precisely known processes are encoded. For these processes either no reliable theoretical prediction 
has been advanced or the different approaches yield different results for the same observable. Therefore, it is necessary 
to verify that these uncertainties do not render our conclusion premature. There are 3 identified uncertainties: the 
<7/VA— >if + cross section, the density dependence of the K + N potential and the lifetime of A in matter if produced 
in a collisions with a sharp energy of two scattering partners. We discuss now how these uncertainties influence our 
results: 

Figure 5, top, shows as well the influence of the unknown NA — > K + AN cross section on this ratio. We confront 
the standard IQMD option (with cross sections for AN interactions from Tsushima et al. [23]) with another option, 
cr(NA) = 3/ '4<t(N N) [26], which is based on isospin arguments and has been frequently employed. Both cross sections 
differ by up to a factor of ten and change significantly the absolute yield of K + in heavy ion reactions but do not 
change the shape of the ratio. 




FIG. 5: Comparison of the measured excitation function of the ratio of the K + multiplicities per mass number A obtained in 
Au+Au and in C+C reactions (Ref. [20]) with various calculations. The use of a hard EoS is denoted by thin (blue) lines, 
a soft EoS by thick (red) lines. The calculated energies are given by the symbols, the lines are drawn to guide the eye. On 
top, two different versions of the NA — > K + AN cross sections are used. One is based on isospin arguments [26], the other is 
determined by a relativistic tree level calculation [27]. The calculation by Fuchs [21] are shown as dotted lines. Middle: IQMD 
calculations with and without KN potential are compared. Bottom: The influence of different options for the life time of A in 
matter is demonstrated. 



The middle part demonstrates the influence of the kaon-nucleon potential which is not precisely known at the 
densities obtained in this reaction. The uncertainties due to the A life time are discussed in the bottom part. Both 
calculations represent the two extreme values for this lifetime [23] which is important because the disintegration of 
the A resonance competes with the K + production. 

Thus we see that these uncertainties do not influence the conclusion that the excitation function of the ratio is quite 
different for a soft EoS potential as compared to a hard one and that the data of the KaoS collaboration are only 
compatible with the soft EoS. The only possibility to change this conclusions is the assumption that the cross sections 
are explicitly density dependent in a way that the increasing density is compensated by a decreasing cross section. It 
would have a strong influence on other observables which are presently well predicted by the IQMD calculations. 

The conclusion that nuclear matter is best described by a soft EoS, is supported by another variable, the dependence 
of the K + yield on the number of participating nucleons ^4 part . The prediction of the IQMD simulations in the 
standard version for this observable is shown in Fig. 6. The top of the figure shows the kaon yield M K + /A palt for 
Au+Au collisions at 1.5 A GeV dbS ct function of the participant number vlp ar t for a soft EoS using different options: 
standard version (soft, KN), calculations without kaon-nucleon interaction (soft, no KN) and with the isospin based 
NA — > NAK + cross section (soft, KN, a*). These calculations are confronted with a standard calculation using the 
hard EoS potential. The scaling of the kaon yield with the participant number can be parameterized by M K + = ^4p art . 

All calculations with a soft EoS show a rather similar value of a - although the yields are very different - while the 
calculation using a hard equation shows a much smaller value. Therefore we can conclude that also the slope value a 
is a rather robust observable. 

The bottom of Fig. 6 shows that a depends smoothly on the compressibility re of the EoS. Whether we include 
the momentum dependence of the nucleon nucleon interaction (with mdi) or not (without mdi) does not change the 
value of a as long as the compressibility is not changed - in stark contrast to the in-plane flow. Again, the measured 
centrality dependence for Au+Au at 1.5 A GeV from the KaoS collaboration [28], a = 1.34 + 0.16, is only compatible 
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with a soft EoS potential. That the dependence of the K + yield on the participants number is a clear signal for the 
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FIG. 6: Dependence of the K + scaling on the nuclear EoS. We present this dependence in form of M K + = j4p art . On the top 
the dependence of M K + /A paTt as a function of A paTt is shown for different options: a "hard" EoS with KN potential (solid 
line), the other three lines show a "soft" EoS, without KN potential and <r(iVA) from Tsushima [27] (dotted line), with KN 
potential and the same parametrization of the cross section (dashed line) and with KN potential and a(NA) = 3/4a(NN). 
On the bottom the fit exponent a is shown as a function of the compressibility for calculations with momentum-dependent 
interactions (mdi) and for static interactions (dashed line) [22]. On the right hand side we compare the energy dependence of 
the fit exponent a for the two EoS. 

EoS also at other beam energies as can be seen from fig. 6 right which displays the excitation function of the fit 
exponent a. Data which follow the curve for a soft equation of state will soon be published [29] 

In conclusion, we have shown That earthbound experiments have now reached a precision which allows to determine 
the hadronic EoS. The experimental results for the two observables which are most sensitive to the hadronic EOS 
are only compatible with theory if the hadronic EoS is soft. This conclusion is robust. Little known input quantities 
do not influence this conclusion. The observation of a neutron star with twice the solar mass seems to contradict 
this conclusion. It points toward a hard hadronic EoS. Both results are quite new and one has not to forget that 
we are comparing non equilibrium heavy ion reactions where about the same number of protons and neutrons are 
present and where mesons and baryon resonances are produced with cold neutron matter in equilibrium. To solve 
this contradiction is certainly a big challenge for both communities in the near future. 
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The stiffness of the hadronic equation of state has been extracted from the production rate of 
K+ mesons in heavy ion collisions around 1 A GeV incident energy. The data are best described 
with a compressibility coefficient k around 200 MeV, a value which is usually called "soft". This 
is concluded from a detailed comparison of the results of transport theories with the experimental 
data using two different procedures: (i) the energy dependence of the ratio of K + from Au+Au and 
C+C collisions and (ii) the centrality dependence of the K + multiplicities. It is demonstrated that 
input quantities of these transport theories which are not precisely known, like the kaon-nucleon 
potential, the AN — > NK + A cross section or the life time of the A in matter do not modify this 
conclusion. 
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Since many years one of the most important challenges 
in nuclear physics is to determine E / A(p,T), the en- 
crgy/nucleon in nuclear matter in thermal equilibrium 
as a function of the density p and the temperature T. 
Only at equilibrium density, po, the energy per nucleon 
E/A(p — po,T = 0) = —16 MeV is known by extrapo- 
lation of the Weizsacker mass formula to infinite matter. 
This quest has been dubbed "search for the nuclear equa- 
tion of state (EoS)" . 

Modelling of neutron stars or supernovae have not yet 
constrained the nuclear equation of state Therefore, 
the most promising approach to extract E/A(p,T) are 
heavy ion reactions in which the density of the colliding 
nuclei changes significantly. Three principal experimen- 
tal observables have been suggested in the course of this 
quest which carry - according to theoretical calculations 
- information on the nuclear EoS: (i) the strength distri- 
bution of giant isoscalar monopole resonances @, 0] , (ii) 
the in-plane sidewards flow of nucleons in semi-central 
heavy ion reaction at energies between 100 A MeV and 
400 A MeV 0] and (iii) the production of K+ mesons m 
heavy ion reactions at energies around 1 A GeV Q. Al- 
though theory has predicted these effects qualitatively, a 
quantitative approach is confronted with two challenges: 
a) The nucleus is finite and surface effects are not neg- 
ligible, even for the largest nuclei and b) in heavy ion 
reactions the reacting system does not come into equilib- 
rium. Therefore complicated non-equilibrium transport 
theories have to be employed and the conclusion on the 
nuclear equation of state can only be indirect. 

(i) The study of monopole vibrations has been very 
successful, but the variation in density is minute. There- 
fore, giant monopole resonances are sensitive to the en- 
ergy which is necessary to change the density of a cold nu- 
cleus close to the equilibrium point pq. According to the- 
ory the vibration frequency depends directly on the force 
which counteracts to any deviation from the equilibrium 
and therefore to the potential energy. The careful analy- 
sis of the isoscalar monopole strength in non-rclativistic 



Q and relativistic mean field models has recently con- 
verged 3 due to a new parametrization of the relativis- 
tic potential. These calculations allow now for the deter- 
mination of the compressibility n = 9p 2 d E ^ p p ' T ^ \p=p 
which measures the curvature of E/A(p,T) at the equi- 
librium point. The values found are around k = 240 
MeV and therefore close to what has been dubbed "soft 
equation of state" . 

(ii) If the overlap zone of projectile and target becomes 
considerably compressed in semi-central heavy-ion colli- 
sions, an in-plane flow is created due to the transverse 
pressure on the baryons outside of the interaction region 
with this flow being proportional to the transverse pres- 
sure. In order to obtain a noticeable compression, the 
beam energy has to be large as compared to the Fermi 
energy of the nucleons inside the nuclei and hence a beam 
energy of at least 100 A MeV is necessary. Compression 
goes along with excitation and therefore the comprcs- 
sional energy of excited nuclear matter is encoded in the 
in-plane flow. It has recently been demonstrated [|| that 
transport theories do not agree quantitatively yet and 
therefore former conclusions [7} have to be considered as 
premature. 

(iii) The third method is most promising for the study 
of nuclear matter at high densities and is subject of this 
Letter. K + mesons produced far below the NN thresh- 
old cannot be created in first-chance collisions between 
projectile and target nucleons. They do not provide suf- 
ficient energy even if one includes the Fermi motion. The 
effective energy for the production of a K + meson in the 
NN center of mass system is 671 MeV as in addition 
to the mass of the kaon a nucleon has to be converted 
into a A to conserve strangeness. Before nucleons can 
create a K + at these subthreshold energies, they have to 
accumulate energy. The most effective way to do this is 
the conversion of a nucleon into a A and to produce in 
a subsequent collision a K + meson via AiV — * NK + A. 
Two effects link the yield of produced K + with the den- 
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sity reached in the collision and the stiffness of the EoS. 
If less energy is needed to compress matter (i) more 
energy is available for the K + production and (ii) the 
density which can be reached in these reactions will be 
higher. Higher density means a smaller mean free path 
and therefore the A will interact more often increasing 
the probability to produce a K + and hence, it has a 
lower chance to decay before it interacts. Consequently 
the K + yield depends on the comprcssional energy. At 
beam energies around 1 A GeV matter becomes highly 
excited and mesons are formed. Therefore this process 
tests highly excited hadronic matter. At beam energies 
> 2 A GeV first-chance collisions dominate and this sen- 
sitivity is lost. 

In this Letter we would like to report that for the third 
approach different transport theories have converged. 
Two independent experimental observables, the ratio of 
the excitation functions of the K + production for Au+Au 
and for C+C [12, El, and a new observable, the depen- 
dence on the number of participants of the K + yield show 
that nucleons interact with a potential which corresponds 
to a compressibility of n < 200 McV in infinite matter in 
thermal equilibrium. This value extracted for hadronic 
matter at densities around 2.5 times the normal nuclear 
matter density is very similar to that extracted at normal 
nuclear matter density. A key point of this paper is to 
demonstrate that the different implementation of yet un- 
solved physical questions, like the NA — > K + AN cross 
section, the KN interaction as well as the life time of the 
nuclear resonances in the hadronic environment do not 
affect this conclusion. 

In order to determine the energy which is necessary 
to compress infinite nuclear matter in thermal equilib- 
rium by heavy ion reactions in which no equilibrium is 
obtained one chooses the following strategy: The trans- 
port theory calculates the time evolution of the quantal 
particles described by Gaussian wave functions. The time 
evolution is given by a variational principle and the equa- 
tions one obtains for this choice of the wave function are 
identical to the classical Hamilton equations where the 
classical two-body potential is replaced by the expecta- 
tion value of the real part of the Bruckner G-matrix. For 
this potential the potential energy in infinite nuclear mat- 
ter is calculated. To determine the nuclear equation of 
state we average this (momentum-dependent) two-body 
potential over the momentum distribution of a given tem- 
perature T and add to it the kinetic energy. Expressed 
as a function of the density we obtain the desired nuclear 
equation of state E/A(p, T). Our two-body potential has 
five parameters which are fixed by the binding energy of 
infinite nuclear matter at po, the compressibility k and 
the optical potential which has been measured in pA re- 
actions. 

Once the parameters are fixed we use the two-body 
potential with these parameters in the transport calcula- 
tion. There is an infinite number of two-body potentials 
which give the same equation of state because the range 
of the potential does not play a role in infinite matter. 



The nuclear surface measured in electron scattering on 
nuclei fixes the range, however, quite well. The uncer- 
tainty which remains is of little relevance here (in con- 
tradiction to the calculation of the in-plane flow which is 
very sensitive to the exact surface properties of the nuclei 
and hence to the range of the potential). 

We employ the Isospin Quantum Molecular Dynam- 
ics (IQMD) [9| approach with the following equations of 
motion: 



and f.-m. 

Pi ~ dn and n ~ dp, 



(1) 



where the expectation value of the total Hamiltonian 
reads as (H) = (T) + {V) with 

i T ) = yiL 

X 1 ^ 2m, 

< y > = T,/Zj 'fi(r,P,t)V^ <fj{r ',p',t)drdr'dpdp\2) 

i j>i 

and fi being the Gaussian Wigner density of nucleon 
i. The baryon-potential consists of the real part of 
the G-Matrix which is supplemented by the Coulomb 
interaction between the charged particles. The former 
can be further subdivided in a part containing the con- 
tact Skyrme-type interaction only, a contribution due 
to a finite range Yukawa-potential, and a momentum- 
dependent part with 



yv = 



VZ 3 

Skyrmc 



^Yuk 



yV 
'Coul 



= tiS(xi - Xj) + t 2 5(xi - x j )p 1 1 (x i ) 
exp{-|f, - Xj\/p} 



Z t Zje 2 



Xj\/p 



f 4 ln 2 (l + t 5 (pi - p j ) 2 )6(x l - Xj) 



(3) 



with Zi, Zj the charges of the baryons i and j. For more 
details we refer to Ref. @. 

We include in this calculation all inelastic cross sections 
which are relevant for the K + production. For details 
of these cross sections we refer to [l(| • Unless specified 
differently, the change of the K + mass due to the kaon- 
nucleon (KN) interaction according to m K (p) = m^(l — 
0.075—) is taken into account, in agreement with recent 

Pa* , to 

self-consistent calculations of the spectral function of the 
K + The A potential is 2/3 of the nucleon potential, 
assuming that the s quark is inert. 

In order to minimize the experimental systematical er- 
rors and the consequences of theoretical uncertainties it is 
better to compare ratios of cross sections rather than the 
absolute values [l2||. We have made sure that the stan- 
dard version of IQMD reproduces the excitation function 
for Au+Au as well as for C+C quite well These ra- 
tios are quite sensitive to the nuclear potentials because 
the compression obtained in the Au+Au collisions is con- 
siderable (up to 3po) and depends on the nuclear equa- 
tion of state whereas in C+C collisions the compression 
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is small and almost independent on the stiffness of the 
EoS. 

Figure [T] shows the comparison of the measured ratio 
of the K + multiplicities obtained in Au+Au and C+C 
reactions [l^l together with transport model calculations 
as a function of the beam energy. 

We see clearly that the form of the yield ratio depends 
on the potential parameters (hard EoS: k = 380 MeV, 
thin lines and solid symbols, soft EoS: k — 200 MeV , 
thick lines and open symbols) in a quite sensible way 
and that the prediction in the standard version of the 
simulation (squares) for a soft and a hard EoS potential 
differ much more than the experimental uncertainties. 
The calculation of Fuchs et al. [H| given in the same 
graph, agrees well with our findings. 

This observation is, however, not sufficient to deter- 
mine the potential parameters uniquely because in these 
transport theories several not precisely known processes 
are encoded. Therefore, it is necessary to verify that 
these uncertainties do not render this conclusion prema- 
ture. Figure [TJ top, shows as well the influence of the 
unknown NA — + K + AN cross section on this ratio. We 
confront the standard IQMD option (with cross sections 
for AN interactions from Tsushima et al. [l(|) with an- 
other option, er(AA) = 3/4a(NN) which is based 
on isospin arguments and has been frequently employed. 
Both cross sections differ by up to a factor of ten and 
change significantly the absolute yield of K + in heavy 
ion reactions but do not change the shape of the ratio. 

The middle part demonstrates the influence of the 
kaon-nucleon potential which is not precisely known at 
the densities obtained in this reaction. The uncertain- 
ties due to the A life time are discussed in the bottom 
part. Both calculations represent the two extreme val- 
ues for this lifetime [Toj which is important because the 
disintegration of the A resonance competes with the K + 
production. 

Thus we see that these uncertainties do not influence 
the conclusion that the excitation function of the ratio is 
quite different for a soft EoS potential as compared to a 
hard one and that the data of the KaoS collaboration are 
only compatible with the soft EoS potential. The only 
possibility to change this conclusions is the assumption 
that the cross sections are explicitly density dependent 
in a way that the increasing density is compensated by a 
decreasing cross section. It would have a strong influence 
on other observables which are presently well predicted 
by the IQMD calculations. 

We would like to add that the smoothness of the excita- 
tion function also demonstrates that there are no density 
isomers in the density regions which are obtained in these 
reactions because the K + excitation function would be 
very sensitive to such an isomeric state [HI]. 

The conclusion that nuclear matter is best described by 
a soft EoS, is supported by another variable, the depen- 
dence of the K + yield on the number of participating nu- 
cleons Ap ar t ■ The prediction of the IQMD simulations in 
the standard version for this observable is shown in Fig. [5] 



+ 



3 



5* 



_Soft 
n i 

-An. 
"Hard* 


i ' i 

t d NA (Tsushima) _ 

d N4 =3/4d NN . 
RQMD, Fuchs et al. - 


H 

- >>> 

"Hard 
h 


f 1 1 h 

• KaoS data _ 

' Soft 
. \ , 

with KN pot 
without KN pot 


£^ 

— N I 

Hard' 


: 1 1 h 

' Soft 

" A Kitazoe 

" 4 phase shift 

i , i 



1.0 1.5 

E lab (GeV) 



FIG. 1: Comparison of the measured excitation function of 
the ratio of the K + multiplicities per mass number A obtained 
in Au+Au and in C+C reactions (Ref. [IH) with various cal- 
culations. The use of a hard EoS is denoted by thin (blue) 
lines, a soft EoS by thick (red) lines. The calculated energies 
are given by the symbols, the lines are drawn to guide the eye. 
On top, two different versions of the NA — > K + AN cross sec- 
tions are used. One is based on isospin arguments [l5|], the 
other is determined by a relativistic tree level calculation [TtJ • 
The calculation by Fuchs [l4| are shown as dotted lines. Mid- 
dle: IQMD calculations with and without KN potential are 
compared. Bottom: The influence of different options for the 
life time of A in matter is demonstrated. 



The top of the figure shows the kaon yield M K + /A part 
for Au+Au collisions at 1.5 A GeV as a function of the 
participant number A part for a soft EoS using different 
options: standard version (soft, KN), calculations with- 
out kaon-nucleon interaction (soft, no KN) and with the 
isospin based NA — > AAA' + cross section (soft, KN, 
a*). These calculations arc confronted with a standard 
calculation using the hard EoS potential. The scaling of 
the kaon yield with the participant number can be pa- 
rameterized by M K + = Ap art . 

All calculations with a soft EoS show a rather simi- 
lar value of a - although the yields are very different - 
while the calculation using a hard equation shows a much 
smaller value. Therefore we can conclude that also the 
slope value a is a rather robust observable. 

The bottom of Fig. [2] shows that a depends smoothly 
on the compressibility k of the EoS. Whether we in- 
clude the momentum dependence of the nucleon nucleon 
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FIG. 2: Dependence of the K + scaling on the nuclear equation 
of state. We present this dependence in form of M K + = Ap art . 
On the top the dependence of M K + /Ap art as a function of 
^Ipart is shown for different options: a "hard" EoS with KN 
potential (solid line), the other three lines show a "soft" EoS, 
without KN potential and a(NA) from Tsushima [l7| (dot- 
ted line), with KN potential and the same parametrization 
of the cross section (dashed line) and with KN potential and 
a(NA) = 3/4a(NN). On the bottom the fit exponent a 
is shown as a function of the compressibility for calculations 
with momentum-dependent interactions (mdi) and for static 
interactions (t^ = 0, dashed line). 



interaction (with mdi) or not (without mdi) does not 
change the value of a as long as the compressibility is not 
changed - in stark contrast to the in-plane flow. Again, 
the measured centrality dependence for Au+Au at 1.5 A 
GeV from the KaoS collaboration pj, a = 1.34 ± 0.16, 
is only compatible with a soft EoS potential. 

This finding is also supported by a more recent anal- 
ysis [l9|, of the in-plane flow which supersedes the 
former conclusion that the EoS is hard [2l[ (made before 
the momentum-dependent interaction has been included 
in the calculations). Due to the strong dependence of 
the in-plane flow on the potential range parameter and 
its dependence on the particles observed these conclu- 
sions are much less firm presently. Comparisons of the 
out-of-plane squeeze of baryons also show a preference 
for a soft equation of state with momentum dependent 
interactions [22j. 

In conclusion, we have shown that the two experimen- 
tal observables which are most sensitive to the potential 
parameters of the nucleon-nucleon interaction are only 
compatible with those parameters which lead in nuclear 
matter to a soft hadronic EoS. This conclusion is robust. 
Uncertainties of the input in these calculations, like the 
KN potential at high densities, the lifetime of the A in 
matter and the A7V — > NK + A cross section do not influ- 
ence this conclusion. The potential parameter k is even 
smaller than that extracted from the giant monopole vi- 
brations. Thus the energy which is needed to compress 
hadronic matter of k < 200 MeV is close to the lower 
bound of the interval which has been discussed in the 
past. 

We would like to thank all members of the KaoS Col- 
laboration for fruitful discussions especially A. Forster, 
P. Senger, C. Sturm, and F. Uhlig. 
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The stiffness of the hadronic equation of state has been extracted from the production rate of 
K + mesons in heavy ion collisions around 1 A GeV incident energy. The data are best described 
with a compressibility coefficient k around 200 MeV, a value which is usually called "soft". This 
is concluded from a detailed comparison of the results of transport theories with the experimental 
data using two different procedures: (i) the energy dependence of the ratio of K + from Au+Au and 
C+C collisions and (ii) the centrality dependence of the K + multiplicities. It is demonstrated that 
input quantities of these transport theories which are not precisely known, like the kaon-nucleon 
potential, the AN — > NK + A cross section or the life time of the A in matter do not modify this 
conclusion. This result from earthbound experiments seems to contradict the recent observation of 
a neutron star with twice the solar mass. 

PACS numbers: 25.75.Dw 

How much energy is needed to compress nuclear matter? The answer to this question, the determination of 
E/A(p,T), the energy/nucleon in nuclear matter in thermal equilibrium as a function of the density p and the 
temperature T, has been considered since many years as one of the most important challenges in nuclear physics. 
This quest has been dubbed "search for the nuclear equation of state (EoS)" . 

Only at equilibrium density, po, the energy per nucleon E/A(p = pa, T = 0) = —16 MeV is known by extrapolating 
the Weizsacker mass formula to infinite matter. Standard ab initio many body calculations do not allow for a 
determination of E / A(p, T) at energies well above the saturation density because the detailed form of the interaction 
among hadrons is not known. Therefore another strategy has been developed. Theory has identified experimental 
observables in nuclear reaction physics or in astrophysics which are sensitive to E/A(p,T). Unfortunately these 
observables depend as well on other quantities which are either unknown or little known (like cross sections with 
resonance in the entrance channel) or difficult to asses theoretically (like the resonance lifetimes in hot and dense 
matter) . It was hoped that comparing many observables for different systems and different energies with the theoretical 
predictions these unknown or little known quantities can be determined experimentally and that finally the dependence 
of the observables on E/A(p, T) can be isolated. 

In astrophysics the nuclear EoS plays an important role in binary mergers involving black holes and neutron stars 
[l[ , in double pulsars Q , in the mass-radius relation of neutron stars [H, H| and in supernovae explosions 5j . For a 
recent review on these topics we refer to Q. Unfortunately, as in nuclear reaction physics, there are always other little 
known processes or properties which have to be understood before the nuclear EoS dependence can be isolated. We 
discuss here as example of the mass-radius relation of neutron stars. Fig. [T] shows the neutron star masses in units 
of the solar mass for different types of binaries. These masses are concentrated at around 1-1.5 solar masses. Fig. 
[5] shows a theoretical prediction of the mass-radius relation for neutron stars using different EoS. Since the nature 
of the interior of neutron stars is not known (in contradiction to what the name suggests) one may suppose that it 
consists of hadrons or of quarks. But even if it consists of hadrons there are speculations that there is a K~ or a 
7r~ condensate or that there are hyperons in equilibrium with nuclear resonances. The same is true if the interior 
consists of quarks. Little known color-flavor locked quark phases may modify the EoS at densities which are reached 
in the interior of the neutron star. For a detailed discussion of all these phenomena we refer to ref. @. The mass to 
radius relation for some EoS are shown in fig. [2] We see that the observed masses of neutron stars are compatible 
with almost all quark or hadron based EoS as long as the radius is unknown. Radii, however, are very difficult to 
measure. Because similar problems appear also for other observables, up to recently the astrophysical observations of 
neutron stars did not help much to narrow down the uncertainty on the nuclear EoS. 

This situation has changed dramatically in the last year with the observation of a neutron star with a mass of two 
solar masses Q- If this observation is finally confirmed the mass/radius prediction of figE] excludes that the interior 
of a neutron star is made by quarks, even a soft nuclear EoS, which will be defined below, will be excluded. 

In heavy ion reactions three observables have been identified which are - according to theoretical calculations - 
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FIG. 1: Measured and estimated masses of neutron stars in radio binary pulsars and in x-ray accreting binaries. Error bars are 
lcr. Vertical dotted lines show average masses of each group (1.62 M0, 1.34 Mq and 1.56 Mq); dashed vertical lines indicate 
inverse error weighted average masses (1.48 M Q , 1.41 M Q and 1.34 Mq). The figure is taken from ref 0] 

sensitive to E/A(p, T) at densities larger than pq: (i) the strength distribution of giant isoscalar monopole resonances 
(ii) the in-plane sidewards flow of nucleons in semi-central heavy ion reaction at energies between 100 A MeV 
and 400 A MeV jlo| and (iii) the production of K + mesons in heavy ion reactions at energies around 1 A GeV [ill ]. 
Although theory has predicted these effects qualitatively, a quantitative approach is confronted with two challenges: 
a) The nucleus is finite and surface effects are not negligible, even for the largest nuclei and b) in heavy ion reactions 
the reacting system does not come into equilibrium. Therefore complicated non-equilibrium transport theories have 
to be employed and the conclusion on the nuclear EoS can only be indirect. 

(i) The study of monopole vibrations has been very successful, but the variation in density is minute. Therefore, 
giant monopole resonances are sensitive to the energy which is necessary to change the density of a cold nucleus 
close to the equilibrium point po- According to theory the vibration frequency depends directly on the force which 
counteracts to any deviation from the equilibrium and therefore to the potential energy. The careful analysis of the 
isoscalar monopole strength in non-relativistic Q and relativistic mean field models has recently converged @ due 
to a new parametrization of the relativistic potential. These calculations allow now for the determination of the 

compressibility re — 9p 2 d E {^p ,T ^ \p=p which measures the curvature of E/A(p,T) at the equilibrium point. The 
values found are around re = 240 MeV and therefore close to what has been dubbed "soft EoS" . 

(ii) If the overlap zone of projectile and target becomes considerably compressed in semi-central heavy-ion collisions, 
an in-plane flow is created due to the transverse pressure on the baryons outside of the interaction region with this 
flow being proportional to the transverse pressure. In order to obtain a noticeable compression, the beam energy has 
to be large as compared to the Fermi energy of the nucleons inside the nuclei and hence a beam energy of at least 
100 A MeV is necessary. Compression goes along with excitation and therefore the compressional energy of excited 
nuclear matter is encoded in the in-plane flow. It has recently been demonstrated fl2| that transport theories do not 
agree quantitatively yet and therefore former conclusions j!3| have to be considered as premature. 

(iii) The third method is most promising for the study of nuclear matter at high densities and is subject of this talk. 
K + mesons produced far below the TV TV threshold cannot be created in first-chance collisions between projectile and 
target nucleons. They do not provide sufficient energy even if one includes the Fermi motion. The effective energy for 
the production of a K + meson in the NN center of mass system is 671 MeV as in addition to the mass of the kaon a 
nucleon has to be converted into a A to conserve strangeness. Before nucleons can create a K + at these subthreshold 
energies, they have to accumulate energy. The most effective way to do this is the conversion of a nucleon into a A 
and to produce in a subsequent collision a K + meson via AN — ► NK + A. Two effects link the yield of produced K + 
with the density reached in the collision and the stiffness of the EoS. If less energy is needed to compress matter (i) 
more energy is available for the K + production and (ii) the density which can be reached in these reactions will be 




FIG. 2: Mass-radius diagram for neutron stars. Black (green) curves are for normal matter (SQM) EoS [for definitions of the 
labels, see Q]. Regions excluded by general relativity (GR), causality and rotation constraints are indicated. Contours of 
radiation radii i?oo are given by the orange curves. The figure is from [J]. 



higher. Higher density means a smaller mean free path and therefore the A will interact more often increasing the 
probability to produce a K + and hence, it has a lower chance to decay before it interacts. Consequently the K + yield 
depends on the compressional energy. At beam energies around 1 A GeV matter becomes highly excited and mesons 
are formed. Therefore this process tests highly excited hadronic matter. At beam energies > 2 A GeV first-chance 
collisions dominate and this sensitivity is lost. 

Here we discuss that for the third approach different transport theories have converged. Two independent experi- 
mental observables, the ratio of the excitation functions of the K + production for Au+Au and for C+C [3 [H|, and 
a new observable, the dependence on the number of participants of the K + yield show that nucleons interact with a 
potential which corresponds to a compressibility of k < 200 MeV in infinite matter in thermal equilibrium. This value 
extracted for hadronic matter at densities around 2.5 times the normal nuclear matter density is very similar to that 
extracted at normal nuclear matter density. A key point here is to demonstrate that the different implementation of 
yet unsolved physical questions, like the A A — > K + AN cross section, the KN interaction as well as the life time of 
the nuclear resonances in the hadronic environment do not affect this conclusion. 

In order to determine the energy which is necessary to compress infinite nuclear matter in thermal equilibrium by 
heavy ion reactions in which no equilibrium is obtained one chooses the following strategy: The transport theory 
calculates the time evolution of the quantal particles described by Gaussian wave functions. The time evolution is 
given by a variational principle and the equations one obtains for this choice of the wave function are identical to the 
classical Hamilton equations where the classical two-body potential is replaced by the expectation value of the real part 
of the Bruckner G-matrix. For this potential the potential energy in infinite nuclear matter is calculated. To determine 
the nuclear EoS we average this (momentum-dependent) two-body potential over the momentum distribution of a 
given temperature T and add to it the kinetic energy. Expressed as a function of the density we obtain the desired 
nuclear EoS E/A(p, T). Our two-body potential has five parameters which are fixed by the binding energy of infinite 
nuclear matter at po, the compressibility K and the optical potential which has been measured in pA reactions. 

Once the parameters are fixed we use the two-body potential with these parameters in the transport calculation. 



4 



There is an infinite number of two-body potentials which give the same EoS because the range of the potential does 
not play a role in infinite matter. The nuclear surface measured in electron scattering on nuclei fixes the range, 
however, quite well. The uncertainty which remains is of little relevance here (in contradiction to the calculation of 
the in-plane flow which is very sensitive to the exact surface properties of the nuclei and hence to the range of the 
potential). 

We employ the Isospin Quantum Molecular Dynamics (IQMD) with momentum dependent forces. All details of the 
standard version of the program may be found in [l6| . The standard version is supplemented for this calculation with 
all inelastic cross sections which are relevant for the K + production. For details of these cross sections we refer to 
[ijj. Unless specified differently, the change of the K + mass due to the kaon-nucleon (KN) interaction according to 
m K (p) = m,Q (1 — 0.075^-) is taken into account, in agreement with recent self-consistent calculations of the spectral 
function of the K + [l8j]. The A potential is 2/3 of the nucleoli potential, assuming that the s quark is inert. The 
calculations reproduce the experimental data quite well as can be seen in fig. [3] where we compare the experimental 
and theoretical K + spectra for different centrality bins for 1.48 AGcV Au+Au. This figure shows as well the influence 
of the K + N potential which modifies not only the overall multiplicity of K + due to the increase of the in medium 
mass but also the spectral form. In order to minimize the experimental systematical errors and the consequences of 

Au(U8AGeV)+Au 0-3. Ofm Au(U8AGeV)+Au 7.24-9.4fm 




PLab<GeV) PLab<GeV> 



FIG. 3: K + spectra for different centrality bins as compared with ( preliminary experimental data 

theoretical uncertainties it is better to compare ratios of cross sections rather than the absolute values [3]. We have 
made sure that the standard version of IQMD reproduces the excitation function for Au+Au as well as for C+C quite 
well These ratios are quite sensitive to the nuclear potentials because the compression obtained in the Au+Au 
collisions is considerable (up to 3/?o) and depends on the nuclear EoS whereas in C+C collisions the compression is 
small and almost independent on the stiffness of the EoS. 

Figure[l]shows the comparison of the measured ratio of the K + multiplicities obtained in Au+Au and C+C reactions 
together with transport model calculations as a function of the beam energy. 

We see clearly that the form of the yield ratio depends on the potential parameters (hard EoS: k = 380 MeV, thin 
lines and solid symbols, soft EoS: k = 200 MeV , thick lines and open symbols) in a quite sensible way and that the 
prediction in the standard version of the simulation (squares) for a soft and a hard EoS potential differ much more 
than the experimental uncertainties. The calculation of Fuchs et al. fljl ] given in the same graph, agrees well with 
our findings. 

This observation is, however, not sufficient to determine the potential parameters uniquely because in these transport 
theories several not precisely known processes are encoded. Therefore, it is necessary to verify that these uncertainties 
do not render this conclusion premature. Figure [5J top, shows as well the influence of the unknown A^A — > K + AN 
cross section on this ratio. We confront the standard IQMD option (with cross sections for AN interactions from 
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FIG. 4: Comparison of the measured excitation function of the ratio of the K + multiplicities per mass number A obtained in 
Au+Au and in C+C reactions (Ref. [3) with various calculations. The use of a hard EoS is denoted by thin (blue) lines, 
a soft EoS by thick (red) lines. The calculated energies are given by the symbols, the lines are drawn to guide the eye. On 
top, two different versions of the NA — > K + AN cross sections are used. One is based on isospin arguments [2C| . the other is 
determined by a relativistic tree level calculation [2l| . The calculation by Fuchs [III] are shown as dotted lines. Middle: IQMD 
calculations with and without KN potential are compared. Bottom: The influence of different options for the life time of A in 
matter is demonstrated. 



Tsushima et al. [T3|) with another option, a(NA) = 3/4a(NN) [2(|, which is based on isospin arguments and has 
been frequently employed. Both cross sections differ by up to a factor of ten and change significantly the absolute 
yield of K + in heavy ion reactions but do not change the shape of the ratio. 

The middle part demonstrates the influence of the kaon-nucleon potential which is not precisely known at the 
densities obtained in this reaction. The uncertainties due to the A life time are discussed in the bottom part. Both 
calculations represent the two extreme values for this lifetime fl7j which is important because the disintegration of 
the A resonance competes with the K + production. 

Thus we see that these uncertainties do not influence the conclusion that the excitation function of the ratio is quite 
different for a soft EoS potential as compared to a hard one and that the data of the KaoS collaboration are only 
compatible with the soft EoS potential. The only possibility to change this conclusions is the assumption that the 
cross sections are explicitly density dependent in a way that the increasing density is compensated by a decreasing 
cross section. It would have a strong influence on other obscrvablcs which are presently well predicted by the IQMD 
calculations. 

The conclusion that nuclear matter is best described by a soft EoS, is supported by another variable, the dependence 
of the K + yield on the number of participating nuclcons A part . The prediction of the IQMD simulations in the 
standard version for this observable is shown in Fig. The top of the figure shows the kaon yield M K + /A palt for 
Au+Au collisions at 1.5 A GeV as a function of the participant number A palt for a soft EoS using different options: 
standard version (soft, KN), calculations without kaon-nucleon interaction (soft, no KN) and with the isospin based 
iVA — ► NAK + cross section (soft, KN, a*). These calculations are confronted with a standard calculation using the 
hard EoS potential. The scaling of the kaon yield with the participant number can be parameterized by M K + = A? art . 

All calculations with a soft EoS show a rather similar value of a - although the yields are very different - while the 
calculation using a hard equation shows a much smaller value. Therefore we can conclude that also the slope value a 
is a rather robust observable. 
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The bottom of Fig. [5] shows that a depends smoothly on the compressibility k of the EoS. Whether we include 
the momentum dependence of the nucleon nucleon interaction (with mdi) or not (without mdi) does not change the 
value of a as long as the compressibility is not changed - in stark contrast to the in-plane flow. Again, the measured 
centrality dependence for Au+Au at 1.5 A GeV from the KaoS collaboration [22|, a = 1.34 ±0.16, is only compatible 
with a soft EoS potential. That the dependence of the K + yield on the participants number is a clear signal for the 
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FIG. 5: Dependence of the K + scaling on the nuclear EoS. We present this dependence in form of M K + = Ap art . On the top 
the dependence of M K + /Apart as a function of Apart is shown for different options: a "hard" EoS with KN potential (solid 
line), the other three lines show a "soft" EoS, without KN potential and a(NA) from Tsushima (dotted line), with KN 
potential and the same parametrization of the cross section (dashed line) and with KN potential and er(jVA) = 3/4ct(AW). 
On the bottom the fit exponent a is shown as a function of the compressibility for calculations with momentum-dependent 
interactions (mdi) and for static interactions (£4 = 0, dashed line). On the right hand side we compare the energy dependence 
of the fit exponent a for the two EoS. 

EoS also at other beam energies as can be seen from fig. [5] right which displays the excitation function of the fit 
exponent a. 

In conclusion, we have shown that the two experimental observables which are most sensitive to the potential 
parameters of the nucleon-nucleon interaction are only compatible with those parameters which lead in nuclear 
matter to a soft hadronic EoS. This conclusion is robust. Uncertainties of the input in these calculations, like the 
KN potential at high densities, the lifetime of the A in matter and the AJV — > NK + A cross section do not influence 
this conclusion. The potential parameter k is even smaller than that extracted from the giant monopole vibrations. 
Thus the compressibility n of hadronic matter of < 200 MeV is close to the lower bound of the interval which has 
been discussed in the past. How this result is compatible with the observation of a neutron star with twice the solar 
mass is presently studied. 

We would like to thank all members of the KaoS Collaboration for fruitful discussions especially A. Forster, P. Sen- 
ger, C. Sturm, and F. Uhlig. 



[1] Sasa Ratkovic, Madappa Prakash, James M.Lattimer Submitted to ApJ, astro -ph /05 1 2 1 36] 

[2] Ph. Podsiadlowsk i, J. D. M. Dewi, P. Lesaffre, J. C. Miller, W. G. Newton, J. R. Stone Mon.Not.Roy.Astron.Soc. 361 

(2005) 1243-1249, |astro-ph/0506566"1 
[3] J.M. Lattimer, M. Prakash Science Vol. 304 2004 (536-542) 

[4] J.M. Lattimer and M. Prakash, ApJ 550 (2001) 426; A.W. Steiner, M. Prakash and J.M. Lattimer, Phys. Lett. B486 
(2000) 239; M. Alford and S. Reddy, Phys. Rev. D. 67 (2003) 074024. 



7 



[5] H.-Th. Janka, R. Buras, F.S. Kitaura Joyanes, A. Marek, M. Rampp 



Procs. 12th Workshop on Nuclear Astrophysics, Ringberg Castle, March 22-27, 2004, |astro-ph/0405289 
Fridolin Weber, Prog.Part.Nucl.Phys. 54 (2005) 193-288 
David J. Nice, Eric M. Splaver, Ingrid H. Stairs, Oliver Loehmer, Axel Jessner, Michael Kramer, James M. Cordes, 
submitted to ApJ, |astro-ph/0508050| 

D.H. Youngblood, H.L. Clark and Y.-W. Lui, Phys. Rev. Lett 84 (1999) 691. 
J. Piekarewicz, Phys. Rev. C69 (2004) 041301 and references therein. 
H. Stocker and W. Greiner, Phys. Reports 137, 278 (1986) and references therein. 
J. Aichelin and CM. Ko, Phys. Rev. Lett. 55 (1985) 2661. 
A. Andronic et al., Phys. Lett. B612 (2005) 173. 
P. Danielewicz, R. Lacey, W.G. Lynch, Science 298 (2002) 1592. 
C. Sturm et al., (KaoS Collaboration), Phys. Rev. Lett. 86 (2001) 39. 
C. Fuchs et al., Phys. Rev. Lett 86 (2001) 1794. 
C. Hartnack et al., Eur. Phys. J. Al (1998) 151. 



E. Kolomeitsev et al., accepted J. of Phys. G, nucl-th/0412037 

C.L. Korpa and M.F.M. Lutz, submitted to Heavy Ion Physics, [micl-th/0404088| 

C. Hartnack and J. Aichelin, Proc. Int. Workshop XXVIII on Gross prop, of Nucl. and Nucl. Excit., Hirschegg, January 
2000 edt. by M. Buballa, W. Norenberg, B. Schaier and J. Wambach; and to be published in Phys. Rep. 
J. Randrup and CM. Ko, Nucl. Phys. A 343, 519 (1980). 

K. Tsushima et al., Phys. Lett. B337 (1994) 245; Phys. Rev. C 59 (1999) 369. 

A. Forster et al., (KaoS Collaboration), Phys. Rev. Lett. 31 (2003) 152301; J. Phys. G. 30 (2004) 393; A. Forster, 
Ph.D. thesis, Darmstadt University of Technology, 2003. 



Recent astrophysical and accelerator based results on the Hadronic Equation of State 
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In astrophysics as well as in hadron physics progress has recently been made on the determination 
of the hadronic equation of state (EOS) of compressed matter. The results are contradictory, 
however. Simulations of heavy ion reactions are now sufficiently robust to predict the stiffness of 
the (EOS) from (i) the energy dependence of the ratio of K + from Au+Au and C+C collisions 
and (ii) the centrality dependence of the K + multiplicities. The data are best described with a 
compressibility coefficient at normal nuclear matter density k around 200 MeV, a value which is 
usually called "soft" The recent observation of a neutron star with a mass of twice the solar mass 
is only compatible with theoretical predictions if the EOS is stiff. We review the present situation. 

PACS numbers: 25.75.Dw 

How much energy is needed to compress nuclear matter? The answer to this question, the determination of 
E/A(p,T), the energy/nucleon in nuclear matter in thermal equilibrium as a function of the density p and the 
temperature T, has been considered since many years as one of the most important challenges in nuclear physics. 
This quest has been dubbed "search for the nuclear equation of state (EoS)" . 

Only at equilibrium density, po, the energy per nucleon E/A(p = po, T = 0) = —16 MeV is known by extrapolating 
the Weizsacker mass formula to infinite matter. Standard ab initio many body calculations do not allow for a 
determination of E/A(p, T) at energies well above the saturation density because the low density many body expansion 
schema ( Bruckner G- matrix) breaks down and therefore the number of contributing terms is exploding. Therefore 
in nuclear reaction physics another strategy has been developed. Theory has identified experimental observables in 
nuclear reaction physics or in astrophysics which are sensitive to E/A(p, T). Unfortunately these observables depend 
as well on other quantities which are either unknown or little known (like cross sections with resonance in the entrance 
channel) or difficult to asses theoretically (like the resonance lifetimes in hot and dense matter). It was hoped that 
comparing many observables for different systems and different energies with the theoretical predictions these unknown 
or little known quantities can be determined experimentally and that finally the dependence of the observables on 
E/A(p,T) can be isolated. 

In astrophysics the nuclear EoS plays an important role in binary mergers involving black holes and neutron stars 
[l| , in double pulsars Q , in the mass-radius relation of neutron stars [1, 0] and in supernovae explosions Q . For a 
recent review on these topics we refer to Q. Unfortunately, as in nuclear reaction physics, there are always other little 
known processes or properties which have to be understood before the nuclear EoS dependence can be isolated. We 
discuss here as example of the mass-radius relation of neutron stars. Fig. [1] shows the neutron star masses in units 
of the solar mass for different types of binaries. These masses are concentrated at around 1-1.5 solar masses. Fig. 
[5] shows a theoretical prediction of the mass-radius relation for neutron stars using different EoS. Since the nature 
of the interior of neutron stars is not known (in contradiction to what the name suggests) one may suppose that it 
consists of hadrons or of quarks. But even if it consists of hadrons there are speculations that there is a K~ or a 7r~ 
condensate or that there are hyperons in equilibrium with nuclear resonances. The same is true if the interior consists 
of quarks. Little known color- flavor locked quark phases may modify the EoS at densities which are reached in the 
interior of the neutron star. For a detailed discussion of all these phenomena we refer to ref. @. We see that the 
observed masses of neutron stars are compatible with almost all quark or hadron based EoS as long as the radius is 
unknown. Radii, however, are very difficult to measure. Because similar problems appear also for other observables, 
up to recently the astrophysical observations of neutron stars did not help much to narrow down the uncertainty on 
the nuclear EoS. 

This situation has changed dramatically in the last year with the observation of a neutron star with a mass of 
two solar masses 0. If this observation is finally confirmed the mass/radius prediction of figE] excludes that the 
interior of a neutron star is made by quarks 0], even a soft nuclear EoS, which will be defined below, will be excluded. 
This is confirmed by the calculation of Maieron Q which uses a MIT bag model or a color dielectric models EoS 
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FIG. 1: Measured and estimated masses of neutron stars in radio binary pulsars and in x-ray accreting binaries. Error bars are 
lcr. Vertical dotted lines show average masses of each group (1.62 M0, 1.34 Mq and 1.56 Mq); dashed vertical lines indicate 
inverse error weighted average masses (1.48 M Q , 1.41 M Q and 1.34 Mq). The figure is taken from ref 0] 



to describe the quark phase. Baldo @ argue that this conclusion may be premature because it depends too much 
on the equation of state of the quark phase. If one replaces the MIT bag model equation of state by that of the 
Nambu - Jona-Lasinio (NJL) Lagrangian under certain conditions (no color conducting phase) larger masses may be 
obtained. The standard NJL Lagrangian lacks, however, repulsion and in view of the momentum cut-off, necessary 
to regularize the loop integrals, and the coupling constants in the diquark sector, which are not uniquely determined 
by the Fierz transformation, quantitative prediction at hi gh q uark densities are difficult in this approach even if 
qualitative agreement with pQCD calculation can be found [10( . 

Simulations of heavy ion reactions have shown that there are three possible observables which are sensitive to 
E/A(p,T) at densities larger than p : (i) the strength distribution of giant isoscalar monopole resonances [TTI. [T^|. 
(ii) the in-plane sidewards flow of nuclcons in semi-central heavy ion reactions at energies between 100 A MeV and 



400 A MeV [13| and (iii) the production of K + mesons in heavy ion reactions at energies around 1 A GeV [14[ . For 
the present status of these approaches we refer to [l5[ . 

Monopole resonances test the nuclear EoS at densities only slightly larger than the normal nuclear matter density. 
Therefore they are of little help if one compares the EoS determined from astrophysics with that extracted from 
nuclear reaction physics. For the in-plane flow the conclusions are not conclusive yet. This is due to the difficulties 
to determine the EoS in heavy ion collisions. An EoS is defined in a thermally equilibrated system but in heavy 
ion collisions equilibrium is not obtained as the momentum distribution of hadrons shows. In addition, nuclei are 
finite size systems where the surface plays an important role. This can easily be seen inspecting the Weizsacker mass 
formula which gives for infinite matter almost twice the binding energy/per nucleon as for finite nuclei. Therefore 
complicated non-equilibrium transport theories have to be employed and the conclusion on the nuclear EoS can only 
be indirect, in determining the EoS for those potentials which give best agreement with the heavy ion results. 

In order to determine the energy which is necessary to compress infinite nuclear matter in thermal equilibrium by 
heavy ion reactions in which no equilibrium is obtained one chooses the following strategy: The transport theory 
calculates the time evolution of the quantal particles described by Gaussian wave functions. The time evolution 
is given by a variational principle and the equations one obtains for this choice of the wave function are identical 
to the classical Hamilton equations where the classical two-body potential is replaced by the expectation value of 
a Skyrme potential. The Skyrme potential is a simple approximation to the real part of the Bruckner G-matrix 
which is too complicated for performing simulations of heavy ion collisions. For this potential the potential energy in 
infinite nuclear matter is calculated. To determine the nuclear EoS we average this (momentum-dependent) two-body 
potential over the momentum distribution of a given temperature T and add to it the kinetic energy. Expressed as a 
function of the density we obtain the desired nuclear EoS E/A(p, T). The potential which we use has five parameters. 




FIG. 2: Mass-radius diagram for neutron stars. Black (green) curves are for normal matter (SQM) EoS [for definitions of the 
labels, see Q]. Regions excluded by general relativity (GR), causality and rotation constraints are indicated. Contours of 
radiation radii i?oo are given by the orange curves. The figure is from [4(. 



Four of them are fixed by the binding energy per nucleon in infinite nuclear matter at po and the optical potential 
which has been measured in pA reactions [161 ] . The only parameter which has been not determined by experiments 
yet is the compressibility k at po- For k < 250 MeV one calls the EoS soft, whereas an EoS is called hard for k > 350 
MeV. Once the parameters are fixed we use the two-body potential with these parameters in the transport calculation. 
There is an infinite number of two-body potentials which give the same EoS because the range of the potential does 
not play a role in infinite matter. The nuclear surface measured in electron scattering on nuclei fixes the range, 
however, quite well. 

The different transport theories give quite comparable results for the bulk part but it is difficult to model the 
surface. (In these simulations there is no surface in the strict sense. Each nucleon contributes to the density by its 
Gaussian wave function and the positions of the hadrons in the course of the reaction determine the surface as well 
as the density gradients.) 

The in-plane flow is caused by the density gradient and hence the numerical value depends on how good the surface 
of the nucleus can be modeled during the reaction. Already small density fluctuations, which arc difficult to control, 
change the value of the in-plane flow considerably. Therefore the second approach, the determination if the EoS by 
measuring the in-plane flow, has not produced conclusive results yet [13] ■ 

The third approach, to measure the EoS by means of the K + yield, depends on bulk properties of matter and 
surface fluctuations have no influence. Here the different transport theories have converged. This was possible due to 
a special workshop at the ECT* in Trento/Italy where the authors of the different simulation codes have discussed 
their approaches in detail and have unified most of the input quantities. The results of this common effort have been 
published in [l8| . As an example we display here the K + p t spectra at midrapidity obtained in the different transport 
theories at different energies. Because with each K + N rescattering collision the slope of the K + spectra changes the 
slope of the pt spectra encodes not only the K + momentum distribution at the time point of production but also 
the distribution of the number of rescatterings. It is therefore all but trivial. Without the KN potential the slopes 
are almost identical and even the absolute yield, which depends on a correct modeling of the Fermi motion of the 
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FIG. 3: Final K + transverse momentum distribution at b—1 fm, |j/ cm < 0.5 1 and with an enforced A lifetime of 1/120 MeV 
(top row without, bottom row with KN potential) in the different approaches [l8l |. 

nucleons, is very similar. If we include the KN interaction which is not identical in the different approaches (see [H|) 
we still observe a very similar slope for most of the programs. 

Due to this progress the simulation programs can now be used to extract up to now theoretically inaccessible 
information like the hadronic EOS [Hj]. Three independent experimental observablcs, the ratio of the excitation 
functions of the K + production for Au+Au and for C+C [lH, l2Cj . the dependence of the K + yield on the number 
of participants and the excitation function of this dependence can be simultaneously reproduced if in these transport 
theories the nucleons interact with potential which yield in infinite matter in equilibrium a compressibility of of the 
EoS of k « 200 MeV. Large compressibility moduli yield results which disagree with all three observablcs. 

This value of k extracted from the K + production which is sensitive to nuclear matter around 2.5po is very similar 
to that extracted by the study of monopole vibrations at p [ll], [l2j . 

It is not sufficient to determine the compressibility modulus. One has to demonstrate as well that its numerical 
value is robust, i.e. that the different implementations of yet unsolved physical questions, like the NA — > K + AN 
cross section, the KN interaction as well as the life time of the nuclear resonances in the hadronic environment do 
not affect its value. 

We employ the Isospin Quantum Molecular Dynamics (IQMD) with momentum dependent forces. All details of the 
standard version of the program may be found in . In addition we have implemented for this calculation all cross 
sections which yield to the production of K + as well as the elastic and the charge exchange KN — * KN reactions. 
The parametrization of the cross section may be found in [l8j]. In the standard version the K + N potential leads to an 
increase of the K + mass in matter, m K (p) = (1 — 0.075^-), in agreement with recent self-consistent calculations 
of the spectral function of the K + [2l[ . The A potential is 2 /3 of the nucleon potential, assuming that the s quark 
is inert. The calculations reproduce the experimental data quite well as can be seen in fig. [5] where we compare the 
experimental and theoretical K + spectra for different centrality bins and for 1.48 AGeV Au+Au. This figure shows 
as well the influence of the K + N potential which modifies not only the overall multiplicity of K + due to the increase 
of the in medium mass but also the spectral form confirming the complexity of the transverse momentum spectrum. 

In order to minimize the experimental systematical errors and the consequences of theoretical uncertainties the 
KaoS collaboration has proposed to study not directly the excitation function of the K + yield but that of the yield 
ratio of heavy to light systems (l9l | . Calculations have shown that ratios are much less sensitive to little known input 
parameters because these affect both systems in a rather similar way. We have shown in fig. 2] that the absolute 
yields are well reproduced in our simulations. Therefore we can use this ratio directly for a quantitative comparison 
with data. The ratio of the K + yields obtained in C+C and Au+Au collisions is quite sensitive to the EoS because 
in Au+Au collisions densities up to 3 po (depending on the EoS) are reached whereas in C+C collisions compression 
is practically absent due to less stopping. 

Figure[5]shows the comparison of the measured ratio of the K + multiplicities obtained in Au+Au and C+C reactions 
[l9t together with transport model calculations as a function of the beam energy [l5| . We see, first of all in the top 




FIG. 4: K + spectra for different centrality bins as compared with (preliminary] experimental data from the KaoS collaboration 



row, that the excitation function of the yield ratio depends on the potential parameters (hard EoS: k = 380 MeV, 
thin lines and solid symbols, soft EoS: K = 200 MeV , thick lines and open symbols) in a quite sensible way and - 
even more essential - that the prediction in the standard version of the simulation (squares) for a soft and a hard EoS 
potential differ much more than the experimental uncertainties. The calculation of Fuchs et al. [2fj| given in the same 
graph, agrees well with our findings. 

This observation is, as said, not sufficient to determine the potential parameters uniquely because in these transport 
theories several not precisely known processes are encoded. For these processes either no reliable theoretical prediction 
has been advanced or the different approaches yield different results for the same observable. Therefore, it is necessary 
to verify that these uncertainties do not render our conclusion premature. There are 3 identified uncertainties: the 
&na^k+ cross section, the density dependence of the K + N potential and the lifetime of A in matter if produced 
in a collisions with a sharp energy of two scattering partners. We discuss now how these uncertainties influence our 
results: 

Figure top, shows as well the influence of the unknown NA — > K + AN cross section on this ratio. We confront 
the standard IQMD option (with cross sections for AN interactions from Tsushima et al. (HI) with another option, 
a(NA) = 3/4a(NN) [23|, which is based on isospin arguments and has been frequently employed. Both cross sections 
differ by up to a factor of ten and change significantly the absolute yield of K + in heavy ion reactions but do not 
change the shape of the ratio. 

The middle part demonstrates the influence of the kaon-nucleon potential which is not precisely known at the 
densities obtained in this reaction. The uncertainties due to the A life time are discussed in the bottom part. Both 
calculations represent the two extreme values for this lifetime [l8[ which is important because the disintegration of 
the A resonance competes with the K + production. 

Thus we see that these uncertainties do not influence the conclusion that the excitation function of the ratio is quite 
different for a soft EoS potential as compared to a hard one and that the data of the KaoS collaboration are only 
compatible with the soft EoS. The only possibility to change this conclusions is the assumption that the cross sections 
are explicitly density dependent in a way that the increasing density is compensated by a decreasing cross section. It 
would have a strong influence on other observables which are presently well predicted by the IQMD calculations. 

The compression which can be obtained in heavy ion reactions depends on the impact parameter or, equivalcntly, 
on the experimentally accessible number of participating nucleons. Therefore by varying the impact parameter we can 
test the EoS at different densities. This dependence should be different for different EoS. This is indeed the case for 
the result of the simulations as seen in Fig. [51 top, where we display the kaon yield M K + /M. pa rt for Au+Au collisions 
at 1.5 A GeV as a function of the participant number A part and for different options: standard version (soft, KN), 
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FIG. 5: Comparison of the measured excitation function of the ratio of the K + multiplicities per mass number A obtained in 
Au+Au and in C+C reactions (Ref. [ltj) with various calculations. The use of a hard EoS is denoted by thin (blue) lines, 
a soft EoS by thick (red) lines. The calculated energies are given by the symbols, the lines are drawn to guide the eye. On 
top, two different versions of the NA — > K + AN cross sections are used. One is based on isospin arguments [23( |. the other is 
determined by a relativistic tree level calculation [24| . The calculation by Fuchs [2(| are shown as dotted lines. Middle: IQMD 
calculations with and without KN potential are compared. Bottom: The influence of different options for the life time of A in 
matter is demonstrated. 



calculations without kaon-nuclcon interaction (soft, no KN) and with the isospin based NA — > NAK + cross section 
(soft, KN, a*). A variation of the KN potential as well as of the K + production cross section change the dependence 
of the K + yield on the number of participants, which can be parametrized by the form M K + = A^ aj . t , only little. 
On the contrary, if we apply a hard EoS, the slope value a changes considerable and is outside of the values which 
are compatible with the experimental results, as shown in the middle part of the figure. In this figure wc display as 
well the insensitivity of our result to the momentum dependence of the nucleon nuclcon interaction. As long as the 
compressibility is not changed the results of our calculations arc very similar independent on whether wc have a static 
or a momentum dependent NN potential. Hence the dependence of the K + yield on the number of participants is 
also a robust variable for the determination of the EoS which supports our earlier conclusion that the EoS is soft. 
Another confirmation that only a soft EoS describes the experimental data is the beam energy dependence of the 
fitted exponent a which is displayed in the right part of fig. [SI The data, which follow the curve for a soft equation 
of state, will soon be published |2q |. 

In conclusion, wc have shown that carthbound experiments have now reached a precision which allows to determine 
the hadronic EoS. The experimental results for the three observables which are most sensitive to the hadronic EOS 
are only compatible with theory if the hadronic EoS is soft. This conclusion is robust. Little known input quantities 
do not influence this conclusion. The observation of a neutron star with twice the solar mass seems to contradict this 
conclusion. It points toward a hard hadronic EoS. Both results are quite new and one has not to forget that we are 
comparing non equilibrium heavy ion reactions where about the same number of protons and neutrons are present 
and where mesons and baryon resonances are produced with cold neutron matter in equilibrium. In addition this 
contradiction depends also on the prediction that the observed star mass excludes the formation of quark matter in 
the interior, a consequence of the suggested EoS of quark matter which is still rather speculative. 
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FIG. 6: Dependence of the K + scaling on the nuclear EoS. We present this dependence in form of M K + = v4p art - On the top 
the dependence of M K + /Apart as a function of Apart is shown for different options: a "hard" EoS with KN potential (solid 
line), the other three lines show a "soft" EoS, without KN potential and a(NA) from Tsushima [24| (dotted line), with KN 
potential and the same parametrization of the cross section (dashed line) and with KN potential and <j(NA) — 3/4:a(NN). 
On the bottom the fit exponent a is shown as a function of the compressibility for calculations with momentum-dependent 
interactions (mdi) and for static interactions (dashed line) [l6|]. On the right hand side we compare the energy dependence of 
the fit exponent a for the two EoS. 



To solve this contradiction is certainly a big challenge for both communities in the near future. 
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Recent astrophysical and accelerator based results on the Hadronic Equation of State 



Ch. Hartnack 1 , H. Oeschler 2 and J org Aichclin 1 * 
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4 rue Alfred Kastler, F-44072 Nantes C'edex 03, France and 
2 Institut fur Kernphysik, Darmstadt University of Technology, 64289 Darmstadt, Germany 

In astrophysics as well as in hadron physics progress has recently been made on the determination 
of the hadronic equation of state (EOS) of compressed matter. The results are contradictory, 
however. Simulations of heavy ion reactions are now sufficiently robust to predict the stiffness of 
the (EOS) from (i) the energy dependence of the ratio of K + from Au+Au and C+C collisions 
and (ii) the centrality dependence of the K + multiplicities. The data are best described with a 
compressibility coefficient at normal nuclear matter density k around 200 MeV, a value which is 
usually called "soft" The recent observation of a neutron star with a mass of twice the solar mass 
is only compatible with theoretical predictions if the EOS is stiff. We review the present situation. 

PACS numbers: 25.75.Dw 

How much energy is needed to compress nuclear matter? The answer to this question, the determination of 
E/A(p,T), the energy/nucleon in nuclear matter in thermal equilibrium as a function of the density p and the 
temperature T, has been considered since many years as one of the most important challenges in nuclear physics. 
This quest has been dubbed "search for the nuclear equation of state (EoS)" . 

Only at equilibrium density, po, the energy per nucleon E/A(p = po, T = 0) = —16 MeV is known by extrapolating 
the Weizsacker mass formula to infinite matter. Standard ab initio many body calculations do not allow for a 
determination of E/A(p, T) at energies well above the saturation density because the low density many body expansion 
schema ( Bruckner G- matrix) breaks down and therefore the number of contributing terms is exploding. Therefore 
in nuclear reaction physics another strategy has been developed. Theory has identified experimental observables in 
nuclear reaction physics or in astrophysics which are sensitive to E/A(p, T). Unfortunately these observables depend 
as well on other quantities which are either unknown or little known (like cross sections with resonance in the entrance 
channel) or difficult to asses theoretically (like the resonance lifetimes in hot and dense matter). It was hoped that 
comparing many observables for different systems and different energies with the theoretical predictions these unknown 
or little known quantities can be determined experimentally and that finally the dependence of the observables on 
E/A(p,T) can be isolated. 

In astrophysics the nuclear EoS plays an important role in binary mergers involving black holes and neutron stars 
[l| , in double pulsars Q , in the mass-radius relation of neutron stars [1, 0] and in supernovae explosions Q . For a 
recent review on these topics we refer to Q. Unfortunately, as in nuclear reaction physics, there are always other little 
known processes or properties which have to be understood before the nuclear EoS dependence can be isolated. We 
discuss here as example of the mass-radius relation of neutron stars. Fig. [1] shows the neutron star masses in units 
of the solar mass for different types of binaries. These masses are concentrated at around 1-1.5 solar masses. Fig. 
[5] shows a theoretical prediction of the mass-radius relation for neutron stars using different EoS. Since the nature 
of the interior of neutron stars is not known (in contradiction to what the name suggests) one may suppose that it 
consists of hadrons or of quarks. But even if it consists of hadrons there are speculations that there is a K~ or a 7r~ 
condensate or that there are hyperons in equilibrium with nuclear resonances. The same is true if the interior consists 
of quarks. Little known color- flavor locked quark phases may modify the EoS at densities which are reached in the 
interior of the neutron star. For a detailed discussion of all these phenomena we refer to ref. @. We see that the 
observed masses of neutron stars are compatible with almost all quark or hadron based EoS as long as the radius is 
unknown. Radii, however, are very difficult to measure. Because similar problems appear also for other observables, 
up to recently the astrophysical observations of neutron stars did not help much to narrow down the uncertainty on 
the nuclear EoS. 

This situation has changed dramatically in the last year with the observation of a neutron star with a mass of 
two solar masses 0. If this observation is finally confirmed the mass/radius prediction of figE] excludes that the 
interior of a neutron star is made by quarks 0], even a soft nuclear EoS, which will be defined below, will be excluded. 
This is confirmed by the calculation of Maieron Q which uses a MIT bag model or a color dielectric models EoS 
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FIG. 1: Measured and estimated masses of neutron stars in radio binary pulsars and in x-ray accreting binaries. Error bars are 
lcr. Vertical dotted lines show average masses of each group (1.62 M0, 1.34 Mq and 1.56 Mq); dashed vertical lines indicate 
inverse error weighted average masses (1.48 M Q , 1.41 M Q and 1.34 Mq). The figure is taken from ref 0] 



to describe the quark phase. Baldo @ argue that this conclusion may be premature because it depends too much 
on the equation of state of the quark phase. If one replaces the MIT bag model equation of state by that of the 
Nambu - Jona-Lasinio (NJL) Lagrangian under certain conditions (no color conducting phase) larger masses may be 
obtained. The standard NJL Lagrangian lacks, however, repulsion and in view of the momentum cut-off, necessary 
to regularize the loop integrals, and the coupling constants in the diquark sector, which are not uniquely determined 
by the Fierz transformation, quantitative prediction at hi gh q uark densities are difficult in this approach even if 
qualitative agreement with pQCD calculation can be found [10( . 

Simulations of heavy ion reactions have shown that there are three possible observables which are sensitive to 
E/A(p,T) at densities larger than p : (i) the strength distribution of giant isoscalar monopole resonances [TTI. [T^|. 
(ii) the in-plane sidewards flow of nuclcons in semi-central heavy ion reactions at energies between 100 A MeV and 



400 A MeV [13| and (iii) the production of K + mesons in heavy ion reactions at energies around 1 A GeV [14[ . For 
the present status of these approaches we refer to [l5[ . 

Monopole resonances test the nuclear EoS at densities only slightly larger than the normal nuclear matter density. 
Therefore they are of little help if one compares the EoS determined from astrophysics with that extracted from 
nuclear reaction physics. For the in-plane flow the conclusions are not conclusive yet. This is due to the difficulties 
to determine the EoS in heavy ion collisions. An EoS is defined in a thermally equilibrated system but in heavy 
ion collisions equilibrium is not obtained as the momentum distribution of hadrons shows. In addition, nuclei are 
finite size systems where the surface plays an important role. This can easily be seen inspecting the Weizsacker mass 
formula which gives for infinite matter almost twice the binding energy/per nucleon as for finite nuclei. Therefore 
complicated non-equilibrium transport theories have to be employed and the conclusion on the nuclear EoS can only 
be indirect, in determining the EoS for those potentials which give best agreement with the heavy ion results. 

In order to determine the energy which is necessary to compress infinite nuclear matter in thermal equilibrium by 
heavy ion reactions in which no equilibrium is obtained one chooses the following strategy: The transport theory 
calculates the time evolution of the quantal particles described by Gaussian wave functions. The time evolution 
is given by a variational principle and the equations one obtains for this choice of the wave function are identical 
to the classical Hamilton equations where the classical two-body potential is replaced by the expectation value of 
a Skyrme potential. The Skyrme potential is a simple approximation to the real part of the Bruckner G-matrix 
which is too complicated for performing simulations of heavy ion collisions. For this potential the potential energy in 
infinite nuclear matter is calculated. To determine the nuclear EoS we average this (momentum-dependent) two-body 
potential over the momentum distribution of a given temperature T and add to it the kinetic energy. Expressed as a 
function of the density we obtain the desired nuclear EoS E/A(p, T). The potential which we use has five parameters. 




FIG. 2: Mass-radius diagram for neutron stars. Black (green) curves are for normal matter (SQM) EoS [for definitions of the 
labels, see Q]. Regions excluded by general relativity (GR), causality and rotation constraints are indicated. Contours of 
radiation radii i?oo are given by the orange curves. The figure is from [4(. 



Four of them are fixed by the binding energy per nucleon in infinite nuclear matter at po and the optical potential 
which has been measured in pA reactions [161 ] . The only parameter which has been not determined by experiments 
yet is the compressibility k at po- For k < 250 MeV one calls the EoS soft, whereas an EoS is called hard for k > 350 
MeV. Once the parameters are fixed we use the two-body potential with these parameters in the transport calculation. 
There is an infinite number of two-body potentials which give the same EoS because the range of the potential does 
not play a role in infinite matter. The nuclear surface measured in electron scattering on nuclei fixes the range, 
however, quite well. 

The different transport theories give quite comparable results for the bulk part but it is difficult to model the 
surface. (In these simulations there is no surface in the strict sense. Each nucleon contributes to the density by its 
Gaussian wave function and the positions of the hadrons in the course of the reaction determine the surface as well 
as the density gradients.) 

The in-plane flow is caused by the density gradient and hence the numerical value depends on how good the surface 
of the nucleus can be modeled during the reaction. Already small density fluctuations, which arc difficult to control, 
change the value of the in-plane flow considerably. Therefore the second approach, the determination if the EoS by 
measuring the in-plane flow, has not produced conclusive results yet [13] ■ 

The third approach, to measure the EoS by means of the K + yield, depends on bulk properties of matter and 
surface fluctuations have no influence. Here the different transport theories have converged. This was possible due to 
a special workshop at the ECT* in Trento/Italy where the authors of the different simulation codes have discussed 
their approaches in detail and have unified most of the input quantities. The results of this common effort have been 
published in [l8| . As an example we display here the K + p t spectra at midrapidity obtained in the different transport 
theories at different energies. Because with each K + N rescattering collision the slope of the K + spectra changes the 
slope of the pt spectra encodes not only the K + momentum distribution at the time point of production but also 
the distribution of the number of rescatterings. It is therefore all but trivial. Without the KN potential the slopes 
are almost identical and even the absolute yield, which depends on a correct modeling of the Fermi motion of the 
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FIG. 3: Final K + transverse momentum distribution at b—1 fm, |j/ cm < 0.5 1 and with an enforced A lifetime of 1/120 MeV 
(top row without, bottom row with KN potential) in the different approaches [l8l |. 

nucleons, is very similar. If we include the KN interaction which is not identical in the different approaches (see [H|) 
we still observe a very similar slope for most of the programs. 

Due to this progress the simulation programs can now be used to extract up to now theoretically inaccessible 
information like the hadronic EOS [Hj]. Three independent experimental observablcs, the ratio of the excitation 
functions of the K + production for Au+Au and for C+C [lH, l2Cj . the dependence of the K + yield on the number 
of participants and the excitation function of this dependence can be simultaneously reproduced if in these transport 
theories the nucleons interact with potential which yield in infinite matter in equilibrium a compressibility of of the 
EoS of k « 200 MeV. Large compressibility moduli yield results which disagree with all three observablcs. 

This value of k extracted from the K + production which is sensitive to nuclear matter around 2.5po is very similar 
to that extracted by the study of monopole vibrations at p [ll], [l2j . 

It is not sufficient to determine the compressibility modulus. One has to demonstrate as well that its numerical 
value is robust, i.e. that the different implementations of yet unsolved physical questions, like the NA — > K + AN 
cross section, the KN interaction as well as the life time of the nuclear resonances in the hadronic environment do 
not affect its value. 

We employ the Isospin Quantum Molecular Dynamics (IQMD) with momentum dependent forces. All details of the 
standard version of the program may be found in . In addition we have implemented for this calculation all cross 
sections which yield to the production of K + as well as the elastic and the charge exchange KN — * KN reactions. 
The parametrization of the cross section may be found in [l8j]. In the standard version the K + N potential leads to an 
increase of the K + mass in matter, m K (p) = (1 — 0.075^-), in agreement with recent self-consistent calculations 
of the spectral function of the K + [2l[ . The A potential is 2 /3 of the nucleon potential, assuming that the s quark 
is inert. The calculations reproduce the experimental data quite well as can be seen in fig. [5] where we compare the 
experimental and theoretical K + spectra for different centrality bins and for 1.48 AGeV Au+Au. This figure shows 
as well the influence of the K + N potential which modifies not only the overall multiplicity of K + due to the increase 
of the in medium mass but also the spectral form confirming the complexity of the transverse momentum spectrum. 

In order to minimize the experimental systematical errors and the consequences of theoretical uncertainties the 
KaoS collaboration has proposed to study not directly the excitation function of the K + yield but that of the yield 
ratio of heavy to light systems (l9l | . Calculations have shown that ratios are much less sensitive to little known input 
parameters because these affect both systems in a rather similar way. We have shown in fig. 2] that the absolute 
yields are well reproduced in our simulations. Therefore we can use this ratio directly for a quantitative comparison 
with data. The ratio of the K + yields obtained in C+C and Au+Au collisions is quite sensitive to the EoS because 
in Au+Au collisions densities up to 3 po (depending on the EoS) are reached whereas in C+C collisions compression 
is practically absent due to less stopping. 

Figure[5]shows the comparison of the measured ratio of the K + multiplicities obtained in Au+Au and C+C reactions 
[l9t together with transport model calculations as a function of the beam energy. We see, first of all in the top row, 




FIG. 4: K + spectra for different centrality bins as compared with (preliminary] experimental data from the KaoS collaboration 



that the excitation function of the yield ratio depends on the potential parameters (hard EoS: k = 380 MeV, thin 
lines and solid symbols, soft EoS: K = 200 MeV , thick lines and open symbols) in a quite sensible way and - even 
more essential - that the prediction in the standard version of the simulation (squares) for a soft and a hard EoS 
potential differ much more than the experimental uncertainties. The calculation of Fuchs et al. [20l ] given in the same 
graph, agrees well with our findings. 

This observation is, as said, not sufficient to determine the potential parameters uniquely because in these transport 
theories several not precisely known processes are encoded. For these processes either no reliable theoretical prediction 
has been advanced or the different approaches yield different results for the same observable. Therefore, it is necessary 
to verify that these uncertainties do not render our conclusion premature. There are 3 identified uncertainties: the 
&na^k+ cross section, the density dependence of the K + N potential and the lifetime of A in matter if produced 
in a collisions with a sharp energy of two scattering partners. We discuss now how these uncertainties influence our 
results: 

Figure top, shows as well the influence of the unknown NA — > K + AN cross section on this ratio. We confront 
the standard IQMD option (with cross sections for AiV interactions from Tsushima et al. (HI) with another option, 
a(NA) = 3/4(t(NN) [23|, which is based on isospin arguments and has been frequently employed. Both cross sections 
differ by up to a factor of ten and change significantly the absolute yield of K + in heavy ion reactions but do not 
change the shape of the ratio. 

The middle part demonstrates the influence of the kaon-nucleon potential which is not precisely known at the 
densities obtained in this reaction. The uncertainties due to the A life time are discussed in the bottom part. Both 
calculations represent the two extreme values for this lifetime [l8[ which is important because the disintegration of 
the A resonance competes with the K + production. 

Thus we see that these uncertainties do not influence the conclusion that the excitation function of the ratio is quite 
different for a soft EoS potential as compared to a hard one and that the data of the KaoS collaboration are only 
compatible with the soft EoS. The only possibility to change this conclusions is the assumption that the cross sections 
are explicitly density dependent in a way that the increasing density is compensated by a decreasing cross section. It 
would have a strong influence on other observables which are presently well predicted by the IQMD calculations. 

The compression which can be obtained in heavy ion reactions depends on the impact parameter or, equivalcntly, 
on the experimentally accessible number of participating nucleons. Therefore by varying the impact parameter we can 
test the EoS at different densities. This dependence should be different for different EoS. This is indeed the case for 
the result of the simulations as seen in Fig. [51 top, where we display the kaon yield M K + /M. pa rt for Au+Au collisions 
at 1.5 A GeV as a function of the participant number A part and for different options: standard version (soft, KN), 



6 




FIG. 5: Comparison of the measured excitation function of the ratio of the K + multiplicities per mass number A obtained in 
Au+Au and in C+C reactions (Ref. [ltj) with various calculations. The use of a hard EoS is denoted by thin (blue) lines, 
a soft EoS by thick (red) lines. The calculated energies are given by the symbols, the lines are drawn to guide the eye. On 
top, two different versions of the NA — > K + AN cross sections are used. One is based on isospin arguments [23( |. the other is 
determined by a relativistic tree level calculation [24| . The calculation by Fuchs [2(| are shown as dotted lines. Middle: IQMD 
calculations with and without KN potential are compared. Bottom: The influence of different options for the life time of A in 
matter is demonstrated. 



calculations without kaon-nuclcon interaction (soft, no KN) and with the isospin based NA — > NAK + cross section 
(soft, KN, a*). A variation of the KN potential as well as of the K + production cross section change the dependence 
of the K + yield on the number of participants, which can be parametrized by the form M K + = A^ aj . t , only little. 
On the contrary, if we apply a hard EoS, the slope value a changes considerable and is outside of the values which 
are compatible with the experimental results, as shown in the middle part of the figure. In this figure wc display as 
well the insensitivity of our result to the momentum dependence of the nucleon nuclcon interaction. As long as the 
compressibility is not changed the results of our calculations arc very similar independent on whether wc have a static 
or a momentum dependent NN potential. Hence the dependence of the K + yield on the number of participants is 
also a robust variable for the determination of the EoS which supports our earlier conclusion that the EoS is soft. 
Another confirmation that only a soft EoS describes the experimental data is the beam energy dependence of the 
fitted exponent a which is displayed in the right part of fig. [SI The data, which follow the curve for a soft equation 
of state, will soon be published |2q |. 

In conclusion, wc have shown that carthbound experiments have now reached a precision which allows to determine 
the hadronic EoS. The experimental results for the three observables which are most sensitive to the hadronic EOS 
are only compatible with theory if the hadronic EoS is soft. This conclusion is robust. Little known input quantities 
do not influence this conclusion. The observation of a neutron star with twice the solar mass seems to contradict this 
conclusion. It points toward a hard hadronic EoS. Both results are quite new and one has not to forget that we are 
comparing non equilibrium heavy ion reactions where about the same number of protons and neutrons are present 
and where mesons and baryon resonances are produced with cold neutron matter in equilibrium. In addition this 
contradiction depends also on the prediction that the observed star mass excludes the formation of quark matter in 
the interior, a consequence of the suggested EoS of quark matter which is still rather speculative. 
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FIG. 6: Dependence of the K + scaling on the nuclear EoS. We present this dependence in form of M K + = v4p art - On the top 
the dependence of M K + /Apart as a function of Apart is shown for different options: a "hard" EoS with KN potential (solid 
line), the other three lines show a "soft" EoS, without KN potential and a(NA) from Tsushima [24| (dotted line), with KN 
potential and the same parametrization of the cross section (dashed line) and with KN potential and <j(NA) — 3/4:a(NN). 
On the bottom the fit exponent a is shown as a function of the compressibility for calculations with momentum-dependent 
interactions (mdi) and for static interactions (dashed line) [l6|]. On the right hand side we compare the energy dependence of 
the fit exponent a for the two EoS. 



To solve this contradiction is certainly a big challenge for both communities in the near future. 
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